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We compare two approaches to non-Markovian quantum evohition: one based on the concept of 
divisible maps and the other one based on distinguishability of quantum states. The former concept 
is fully characterized in terms of local generator whereas it is in general not true for the latter one. 
A simple example of random unitary dynamics of a qubit shows the intricate difference between 
those approaches. Moreover, in this case both approaches are fully characterized in terms of local 
decoherence rates. As a byproduct it is shown that entropy might monotonically increase even for 
non-Markovian qubit dynamics. 



PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 

Quantum dynamics is represented by the dynamical 
map, that is, a family of completely positive and trace 
preserving maps At [t > 0) such that Ao = 11. If p 
is an initial state of the system then pt = At(/c) de- 
fines evolution of p. This description provides general- 
ization of unitary evolution At(p) = UtpU} with unitary 
Ut = e~*^*, where H represents the Hamiltonian of the 
system. Any departure from unitary evolution signals 
the nontrivial interaction of the system with an environ- 
ment which is responsible for decoherence and dissipation 
processes [ij in open quantum system. The traditional 
approach to the dynamics of such systems consists in ap- 
plying a suitable Born-Markov approximation leading to 
the celebrated quantum Markov semigroup H-Q] which 
neglects all memory effects. Recent theoretical activity 
and technological progress show the importance of more 
refine approach based on non-Markovian evolution. Non- 
Markovian quantum dynamics becomes in recent years 
very active field of both theoretical and experimental 
research (see e.g. recent papers Q-H^l and references 
therein). 

We shall call quantum evolution Markovian if the cor- 
responding dynamical map At is divisible [l6j . Let us 
recall that At is divisible if At = Vf^sAg and Vt,s is com- 
pletely positive and trace preserving for all t > s, that 
is, it gives rise to 2-parameter family of legitimate prop- 
agators. The essential property of Vf.s is the following 
(inhomogeneous) composition law Vt^s Vs,u = Vf,„ for all 
t > s > u. It should be stressed that Markovian dy- 
namics (divisible map) is entirely characterized by the 
properties of the local in time generators Lt, that is, if 
At satisfies At = LfAt, then At corresponds to Markovian 
dynamics if and only if Lt has the standard form 0, Q 
for alH > 0, that is, 

Ltp = -z[Hit),p]+J2 (va.{t)pV^it) - l{V^it)V^it),p}] , 

with time dependent Hamiltonian H{t) and noise oper- 
ators Va{t). A different approach to Markovianity was 



proposed by Breuer et. al. [T^: authors of [T3l define 
non-Markovian dynamics as a time evolution for the open 
system characterized by a temporary fiow of information 
from the environment back into the system and manifests 
itself as an increase in the distinguishability of pairs of 
evolving quantum states: 
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where \A\\i = TrvA^A denotes the trace norm. Accord- 
ing to 17| the dynamics At is markovian iff <j{pi, P2',t) < 
for all pairs of states pi,p2 and t > 0. Contrary to 
divisible map this approach does not correspond to any 
composition law and is not characterized by the prop- 
erties of local generator. It turns out that condition 
a{pi,p2',t) < is less restrictive than requirement of 
complete positivity for Vt^s and one can construct At 
which is non-Markovian (not divisible) but still gives rise 
to the negative flow of information (see [l3l - [T5| ). 

In this report we provide a simple example of qubit 
dynamics for which the condition a{pi,p2;t) < is fully 
characterized by local generator. Consider the following 
random unitary dynamical map 



Atp 



Q=0 



Pa{t)(TaP<ya 
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where uo = and (Ji,U2t<^3 are Pauli matrices, and 
Pa{t) is time-dependent probability distribution such 
that po(0) = 1 which guarantee that Ao = 1. This map 
was recently investigated in [2^ . It is clear that one may 
consider ([2]) as a family of Pauli channels. To answer the 
question wether At is Markovian one has to analyze the 
corresponding local generator. To find Lt = AtA^"^ one 
has to compute the inverse map A^^. Let us observe that 



At(crQ) = Xa{t)cra 



(3) 
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where the time-dependent eigenvalues are given by 

3 

13=0 

with Hap being the Hadaniard matrix 



H = 



for k = 1, 2, 3 and t > 0. The inverse relations yield 



(4) 
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Note that Ao(t) = 1 and |Afc(t)| < 1 for fc = 1,2,3. It is 
therefore clear that 



where 



Kit) 



Ac.(t) 



(6) 
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and hence in particular i-ia{t) = due to Xo(t) = 1. In- 
troducing 



1 ^ 

7a W = j'^Hapfifjit) 



(8) 



f3=0 

the local generator Lt is given by 
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(9) 



a=0 

Finally, observing that 

3 1 
4 

a=0 a,l3=0 

one arrives at the standard form 
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^tip) = 7fc(i)(o-fcPcrfe - p) 



(10) 



k=l 



Hence 



Proposition 1 The random unitary dynamics (0) is 
Markovian iff 

7iW>0, 72W>0, 73W>0, (11) 
for all t > 0. 

Note, that ([TT]) yields highly nontrivial conditions for 
probability distribution pa (t) : 

E^^-^Sf^^I >0. (12) 



J2a=oHpaPa(t) 



Po{t) = 


1 
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[1 + Ai2(t)- 


f ^13(0 


+ ^23(0] , 


Pl{t) = 


1 

4 


[1-Ai2(t)- 


-^13(0 


+ A23{t)] , 


P2{t) = 


1 

4 


[1-Ai2(t)- 


^A3(t) 


-A23{t)] , 


Pzit) = 


1 
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[1 + ^12(0- 


-^13(0 


-^23(0] , 


where Aij (t) = 


e~ 


2[r.(t)+r,(t)] 


and wc 


introduced 



^k{t) = / n{r)dT . 
Jo 

Let us recall [29 that if Af is Markovian then 



^l|A*(X)||i<0, 
for all Hermitian X. Taking X = ak one obtains 



(13) 



|l|A*(a..)||i = ||Afe(t)|||a,||i<0, (14) 

and hence Markovianity of random unitary evolution ^ 
implies 



^|Afc(i)| <0, 



(15) 



for k = 1,2,3. One easily finds the following relations 
Ai(t) ^ ,-2[r.(t)+r,it)] ^ 



Ht) 
Ht) 

and hence (|15l) implies 



,-2[ri(t)+r3(t)] 



Xsit) - e-2[ri(*)+ri(*)] 



7i(i)+72(t) > 
7i(i)+73(i) > 
72(i)+73(i) > 



(16) 



for all t > 0. We stress that the above conditions are 
only necessary but not sufficient for Markovianity: pT|) 
imply (jl6p but the converse is of course not true. 
Now comes our main result 

Proposition 2 The random unitary dynamics (0i sat- 
isfies a{pi,p2;t) < if and only if conditions il6]) are 
satisfied. 

Proof: let us consider ||At||i, where Af = At (A) and A = 
Pi — p2- Since A is traceless one has A = J2k=i ^kO'k with 
real Xk- Note that A^ = (X^Li^D^^- ^o^' At(A) = 
X]fc=i ^k{t)xkO'k is traceless as well, and hence 

||At(A)||i=Trv/[A;(AF = Tr[e(t)I]=2e(t) , (17) 
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where £_{t) = \>'k{t)\'^xl > 0. Observe that 



illlAII = ^ 
2dt" dt 2£,{t)f^^"'''dV 



It is therefore clear that ^||At||i < for all pairs pi and 
P2 if and only if ^|Afc(t)| < for fc = 1,2,3 which is 
equivalent to ([T6| . 

Consider now a general (not necessarily traceless) Her- 
mitian X = xoI + A with Tr A = 0. Due to At (I) = I one 
has Xt = h.t{X) = xqI + At. Let 6+{t) and 5^{t) denote 
the eigenvalues of At such that 5+{t) ~ —5-{t) > 0. The 
corresponding eigenvalues of Xf read x± (t) = Xo±S+ (t) 
and hence 

\\Xt\\, = \xo + S+{t)\ + \xo-S+{t)\. 

Assuming xo > (for xq < the argument is analogous) 
one finds 



Xq - S+{t) \ d 

'-\xo-5^{t)\idt^^^'^^ 



and hence ^ ||A'f||i = if (5+(t) < xq, and ^ H-'^tlli 
2^ (5+(0 if S+it) > xo. For xo = one has ^ ||At||i 
^ ||Xt||i = 2|^ 5+(i) < which finally implies 



dt 



|At(X)||i<0, 



for arbitrary Hermitian X. 

Interestingly, similar result holds for the entropy S{pt). 
Note that At (I) = I and hence S{pt) > S{p). Now, if At 
is Markovian then -^S{pt) > for all initial states p. 



iSiPt) > 0. 



Proposition 3 // conditions U6\) are satisfied then 
A • 



and hence 4fS{pt) > if i+(t) < 0. One finds 



d , , 



g-[ri(t)+r2(t)] 



Ait) 



where 



A{t) = -[7i(t)+72(t)]u2 

- [7i(0+73(t)]a^e2[r^(*)-r3(*)] 

- [72(0 + 73(0] fe'e^'"^^*^-"^^*^! , 

and hence A{t) < whenever (fT6|) is satisfied. It is clear 
that similar result holds if we replace von Neumann en- 
tropy by Renyi or Tsallis generalized entropies (26| . 
Consider for example 



71 (i) = 72(t) = 1 , 73 (t) = sint 



(21) 



This choice is perfectly legitimate since ri{t) = T2{t) = 
t > Q and T3{t) = 1 — cost > 0. The corresponding evo- 
lution is non-Markovian due to the fact that 73 (t) can 
be negative but the conditions ((T6)) are satisfied. Hence, 
even if one time-dependent decoherence rate is periodic 
and takes strictly negative values both information flow 
and von Neumann entropy behaves perfectly monotoni- 
cally as in the case of Markovian dynamics. 

Interestingly, if there are at most two decoherence 
channels, that is, at most two 7fc(t) are non- vanishing, 
then (HH) imply that 7^ (t) > and hence and 
are equivalent. However, the relation between 7fe(t) and 
probability distribution Pa{t) is still quite complicated. 
If there is only one decoherence channel, say fcth, that is 
7;(t) for ? 7^ fc, then 



and pi{t) = for I ^ k. In this case 



To prove it observe that decomposing 



u w 

w* —u 



with w = a + ib, one finds for the time evolution 



^^ = l]I + ,-[r.(*)+r.(*)][^ n w{t) 
2 \'w*(t) -u 



(18) 



(19) 



where w{t) = e'^^lt) [gri(t) ^ _^ jgr2(t) ^.^j-^) > 

X- [t) be the corresponding eigenvalues of pt 



x±{t) = i ± e-PiW+r=Wl ^u^ + \w{t)\^ 



One has 



(20) 



Pk{t) 



Po{t) 



l-2pfc(t) 2po(0-l 



(22) 



and hence Markovianity is equivalent to pq (t) = —pk (t) < 
0. This way we can fully characterize (non)Markovianity 
of the bit-fiip (fc = 1), phase-flip (k = 2), and phase- 
damping (fc = 3) dynamics. If pi(t) ~ P2{t) = P3(^) — 
i[l-Po(t)], then 



At=p(t)]l2 + [l-p(t)]AD 



(23) 



where p{t) = i[4po(i) — 1] and Ad is a completely de- 
polarizing channel corresponding to Pa{t) = \- Note 
that in this case Wt = {"^-^ Kt)P2 defines a family of 2- 
qubit Werner states satisfying initial condition Wo = P2 
{P2 = 5 J N*)Oj I denotes canonical maximally en- 
tangled state) . Such special qubit channels were recently 
analyzed in Ref. [l^]. One has \k{t) = p(t) for fc = 1, 2, 3 
and hence jkit) = ~'P{t)/[4:p{t)]. It leads to the following 
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local generator Lt{p) = Y.l^ii'^kP crk - p) which 

may be rewritten as follows 

L, = -M[A^_]1,] , (24) 

which shows that the corresponding dynamics is Marko- 
vian iff ^p{t) < 0. 

It is clear that ^ may be generalized for arbitrary 

> 2 as follows 

Up) - E Pkimkipul, (25) 

k,l=Q 

where pki{t) is time-dependent probability distribution, 
and Uki are generalized spin operators defined by 

with Lo = e^'^'/^ (we add modulo N). It is legitimate 
dynamical map provided poo (0) = 1 which guarantee that 
Ao = 1. Using 

UklUrs - Uj'^'Uk+r^l+s , Ul = C^'='C/_fe,_z , (27) 

one finds Kt{Uij) — Xij{t)Uij , with time-dependent 
eigenvalues 

AT-l 

A,,(t) = E Pki{t)u;'''-'' ■ (28) 

fc,/=0 

Note, that Xoo{t) = 1 and \X^J{t)\ < 1 for all t > 0. 
Clearly one has Ay (0) = 1 for all «, j = 0, . . . , N — 1. 
Introducing two NxN matrices X{t) and p{t) with matrix 
elements Ay(t) and Pij(t), respectively, one can rewrite 
([28|) in the following compact form 

A(<) ^ Fp{tfF^ : (29) 

where F denote the matrix of the discrete Fourier trans- 
form, that is, Fii = a;"*'. Using F~^ = j^F^ one finds 
for the inverse relation 

Pitf -^Fnit)F . (30) 

The corresponding time-dependent local generator Lt 
reads 

Lt{p)=J2'^^^'^*HU^,pUl-p) , (31) 

where = E(.j)5:^(o,o) ' ^nd 

l{tf = ^F^Mt)F , (32) 

and pij (<) are eigenvalues of Lt which are related to Xij (t) 
via 

Xij{t) = exp [ / /iy (T)(iT] . 
Jo 

We can summarize 



Proposition 4 T/ie random unitary dynamics I125\) is 
Markovian (At is divisible) if and only if jij{t) > for 
each pair {i,j) 7^ (0, 0) and t > 0. 

In particular taking Pij{i) = [1 — poa{t)]/[N'^ — 1] for 
(i, j) ^ (0,0) one finds At p(i)]ljv + [1 -p(<)]Ad, where 
p{t) = [Af%o(i) - 1]/[A^^ - 1] and Ad is a completely 
depolarizing channel corresponding to Pa{t) ~ ■ In 
this case Wt = (Hat (81 At)P^ defines a family of isotropic 
states in (8) with initial condition Wq = . 
One has Xki{t) = p[t) for (fc, Z) ^ (0,0) and hence 
'-)/j,;(f) = —p(t)/[N'^p{t)]. It leads to the following local 

generator Lt{p) = -j^r^ UkA^i^'P^ki " P] ^^ich may 

be rewritten as Lt = [An — Hat], which shows that 

the corresponding dynamics is Markovian iff ^p{t) < 0. 
For 2 it reduces to ((M)) . 

It is clear [2^ that if At is Markovian then 

||A,,(t)|<0, (33) 

for i,j = 1,...,A^ — 1. Clearly, one obtains a nontrivial 
condition only for {i,j) ^ (0,0). Using X^^^Lo 7mn(0 = 
0, the above formula imply the following conditions for 
the time-dependent decoherence rates 

{t)>0, (34) 

m,n 

where amn;ij = 1 ^ Rca;^"~™. For N = 2 these condi- 
tions reduce to If n = 3 the formula provides 
the following set of conditions 

7io(t) + 711 (0 + 712 (0+72o(t)+ 721 (i) + 722(0 > 0, 
701 (i) + 702 (0 + 711 (0 + 712 (t)+ 721 (i) + 722(0 > 0, 

701 (t) + 702(0 + 710(0 + 712(0 + 720 (0 + 721(0 > 0, 

701 (0 +702 (0 + 710 (0 +711 (0 +720 (0 + 722(0 > 0. 

These are of course much more difficult to handle than 
simple condition (|16p . However, preliminary numerical 
analysis shows that these conditions might be also suffi- 
cient for the negative information flow. 

In conclusion we provide a simple example of qubit ran- 
dom unitary dynamics ^ and show that two concepts of 
Markovian dynamics based on divisible dynamical maps 
and negative fiow of information diverge. Interestingly 
both concepts arc fully characterized in terms of local 
generator: divisibility by and negative flow of infor- 
mation by (Hn]). We stress that it is the first example 
of quantum dynamics which provide full characterization 
of the Breuer et. al. condition [l^l in terms of local in 
time generator Lf. Interestingly, condition ([TS]) is also 
sufficient for the monotonicity of von Neumann entropy 
-3jS{pt) > for all initial states p. Hence, the entropy 
may monotonically increase even for non-Markovian ran- 
dom unitary qubit dynamics. Our discussion can be eas- 
ily generalized for arbitrary N > 2. We conjecture that 
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conditions (p4|) arc sufficient for the negative information ence Center project DEC-2011/03/B/ST2/00136. We 
flow. This point, however, deserves further investigation. thank Ifiigo Luis Egusquiza for his remarks. 
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